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Abstract 

A proposed experiment to measure Kl — ► 7r°z^ at the AGS (E926) is described. 
It adopts a unique but general strategy of rare decay measurements with two inde- 
pendent rejection criteria, which allows reliable measurements of background levels. 
The method employs hermetic veto and full kinematic reconstruction using kaon 
time of flight and kinematic reconstruction of 7r° . Backgrounds are suppressed to a 
level well below an anticipated signal in the range 3 ± 2 x 10~ n . 
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1 Introduction 



1.1 Physics Motivation 

As discussed by the previous speakers, the rare decay Kl — * tx^vv is unique among 
potential SM observables; it is dominated by direct CP violation |]J and is entirely gov- 
erned by short- distance physics involving the top quark. As a consequence of unprece- 
dented theoretical precision and anticipated experimental accessibility, a measurement of 
K° L — > 7T°uu can unambiguously test the SM origin of CP violation and ultimately yield 
the most accurate determination of the CKM CP violating phase 77. This rare decay 
mode therefore provides an exceptional and unique opportunity for making significant 
progress in our understanding of flavor-dynamics and CP violation. It is competitive 
with and complementary to future measurements in the B meson system. Absence of 
K® — ► ir°uu within the expected range of about (3 ± 2) • 10~ n or a conflict with other 
CKM determinations would certainly indicate new physics. 

1.2 A strategy of the rare decay measurement 

In this section, a unique but general strategy for rare decay measurements which was 
found to be essential in observing K + — > 7i + uu in E787, is described. The main challenge 
in such a very rare decay search is to control the systematic error in estimating the level 
of backgrounds. We have to make sure that there are no very rare software/hardware 
problems or unexpected rare background mechanisms that are not taken into account in 
estimating the background. In general, these potential problems are not well represented 
by Monte Carlo simulations, and one has to go back to the data to systematically study 
them. E787 analysis uses a tool called a "bifurcated analysis" along with a method to 
protect against bias due to small statistics. As will become clear, it is essential to have 
two independent rejection criteria ( or cuts ) for each background processes. 

1.2.1 Bifurcated analysis 

The bifurcated analysis of K + — > tt + uu in the E787 experiment consists of two built-in 
independent selection criteria which are used to measure each background. Figure [I] shows 
K + — ► 7r + 7T° background measurement as an example. In this example, the photon veto, 
which rejects ir° — > 77 decays, and the charged particle kinematics (range, energy and 
momentum), which rejects the monochromatic 7r + , are taken as the two independent cuts. 
Because these two cuts address different particles in the final state, they are independent 
in principle. The scatter plot on the left in Figure [l] shows the energy versus range 
distribution when the photon veto cut is loosened. With this sample enhanced by a 
factor of 50, the rejection of the kinematic cut can be measured. The histogram on the 
right shows the momentum distributions before (the largest histogram) and after (the 
smallest histogram) the photon veto. The momentum peak is suppressed by a factor of 
10 6 . It should be pointed out that the momentum line shape did not change after this 
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Figure 1: An example of a bifurcated analysis in the E787 K + — > 7r + 7r° background 
analysis. 



suppression, providing confidence that the two criteria are independent. All the kinematic 
variables used in the analysis are tested this way to see if there are any correlation with 
photon veto left. 

The bifurcated analysis is a powerful way of measuring the background levels and their 
systematic uncertainties due to correlations. Because one can enhance the background 
by a large number, it also allows us to examine rare mechanisms and to develop detailed 
analysis cuts to reject them. Figure |^ shows time differences between the charged track 
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Figure 2: The time difference between the charged track and extra junk energies for the 
K + — > ir + 7T° kinematic peak events. 

and extra energy for the K + — > 7r + 7r° kinematic peak events in the 3 types of photon veto 
systems. The barrel veto is a lead/plastic scintillator sandwich, the range stack is a stack 
of plastic scintillators, and the endcap is made of pure Csl crystals. The solid lines show 
the data and the dashed lines show measured accidental rate from K + — > fi + u events. For 
the barrel veto and range stack, there are non-Gaussian tails at larger times, which are 
understood as neutrons produced by photonuclear reactions (Figure §). Neutrons from 
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photonuclear reactions leave late energy through n-p scattering in the plastic scintillation 
counters in the range stack and barrel veto. Using these data, photon veto cut positions 
in time versus energy space were optimized. 
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Figure 3: Photonuclear reaction cross section as a function of photon energy. The giant 
dipole resonance of nuclei is at around 20MeV, and the delta resonance of the nucleon is 
at around 300MeV. The cross section converges to an asymptotic value at high energy. 



1.2.2 Protection against bias due to small statistics 

From time to time, we see experimental evidence of a several standard deviation (a) effect 
appear and then disappear as more statistics are accumulated. This can be understood 
considering only statistical effects, once we recognize the fact that one measurement out 
of a large number of experiments being performed should show such a large statistical 
fluctuation. The same thing would happen even in a single peak search experiment, such 
as Higgs search: Suppose we look for a peak in a distribution which spans a range 300 
times larger than the width of the expected signal peak, we would expect on average one 
3o" peak, whose probability is 0.27% (300 x 0.27% ~ 1). Confusion is caused when the 
number of trials, e.g. 300 of peak search regions, is not taken into account. This is often 
referred to as analysis bias. 

The same bias shows up in rare decay analyses. Suppose we are left with one candidate 
event and found that one signal selection criterion, which is 99% efficient, removes it. We 
might conclude that the event is a background. However, if 100 such 99% efficient criteria 
are tried before finding the one criterion that removes the event, it is likely we would 
succeed in finding one even for a signal event (100 x 1% = 1). Again, what is important 
is to take the number of trials into account. In reality, we often look at an event display 
and various quantities to make a physics judgment on the event, in which procedure it 
is very difficult to quantify how many trials are performed. One way, which may be the 
only way, to protect against this analysis bias is to avoid such procedures in which one 
cannot quantify the number of trials. 
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There were three steps taken by the rare decay analysis, each of which required pro- 
tection against potential bias: 

• Candidate event selection (Blind analysis) 

If we examine individual candidate events and develop selection cuts, we introduce 
analysis bias as discussed above. The way to avoid this is to develop cuts without 
examining the candidate events ("blind analysis"). The bifurcated analyses provide 
large data samples which are good representations of background. These samples 
allow us to develop analysis cuts without examining candidate events. 

• Background measurement (Training/Test samples) 

If there are no events left (upper limit only) at the end, the blind analysis might be 
sufficient. However, in order to observe the signal, one needs to avoid bias in the 
background measurement, to prove if the remaining events are signal or background. 
Again, when we examine each of the remaining background events, making many 
trials to find effective cuts and adopt only the ones that eliminate the remnant, we 
introduce a bias, namely we underestimate the background level. Because bifurcated 
analysis provides large background samples, we can afford to split them into two 
and use one for developing cuts (training sample) and the other for measuring the 
resulting background levels (test sample) to avoid the potential bias. 

• Detailed examination of the candidates (Signal/background likelihood) 

It is impossible to fully prove that there is nothing missed at all in the analysis 
at the end. In the case of bifurcated analysis, there may be a hidden correlation 
between two "independent" cuts which does not show up in the correlation tests as 
described in the previous section. We would eventually be forced to examine the 
candidate event to see if there is anything wrong with it. As discussed above, we 
could always find a distribution in which the candidate event is near the edge even 
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Figure 4: Signal likelihood function of E787 in the signal box. The shaded histogram 
shows the expected background distribution scaled up by a factor of 10. The candidate 
event is found to be in the cleanest region with a factor of 10 more rejection. 
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for the real signal event, if we look at enough quantities (enough trials). In order 
to examine the candidate event in an unbiased way, a signal likelihood function was 
prepared beforehand in the E787 analysis (Figure f|). The large background sample 
from bifurcated analysis provided enough statistics even inside the signal region. 
There was more confidence in the signal, because the event was in the cleanest 
region in the signal box with a factor of 10 more background rejection. It should be 
pointed out that this last step was feasible only because the experiment had extra 
background rejection capability (redundancy) beyond the final cut position. 



2 Basic concept of the experiment 

Figure |5| shows the ir° energy distribution in the Kl center of mass frame (E* ) for the 
Kl — > tt°uu and Kl — > 7r°7r° (K® 2 ) decays. The main background comes from the CP 
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Figure 5: Energy spectrum of 7r° in the Kl center of mass frame. 

violating decay Kl — > 7r°7r° where 2 of the photons are missed. Even with the branching 
ratio of 9 x 10~ 4 and photon detection inefficiency of 10 -4 per photon (the value that is 
achieved for higher energy photons in E787), the background appears at a single event 
sensitivity of (9 x 10~ 4 ) • (10~ 4 ) 2 -4 C 2 = 5.4 x 10 -11 . This is larger than the central value 
of the Standard Model prediction of 3 x 10~ n . As we have seen in the previous section, it 
is very difficult to measure the background with only a photon veto cut. Thus, to increase 
the probability that the source of observed 'V plus nothing" events is truly the ir°uu mode 
another handle using kinematics is needed. By measuring the decaying K° L momentum 
using time-of-flight, one can transform the 7r° momentum back to the Kl center of mass 
system on an event-by-event basis, and reject the monochromatic K® 2 peak. With these 
two independent cuts of kinematics and photon veto, we can get sufficient background 
rejection and make a reliable measurement of the background. 

Figure || shows a schematic view of the detector. A bunched Kl beam produced by 
the pulsed AGS proton beam decays in the decay region which is hermetically surrounded 
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Figure 6: A schematic view of the detector. 



by photon detectors. The directions of the photons are measured by a preradiator and 
the total energy is measured by an electromagnetic calorimeter. In the following sections, 
more detailed descriptions of the kinematic criteria and photon veto are given. 

2.1 Full kinematic reconstruction and particle identification 

Since the 3-body 7r° spectrum is the only observable in K\ — > 7r°z/P, the most effec- 
tive strategy for a clear measurement is to fully reconstruct the 7r° in the Kl center of 
mass frame. This can only be efficiently and unambiguously accomplished by completely 
measuring the kinematics of the decay photons including time, position, angle and en- 
ergy and by determining the Kl momentum by time-of- flight from the production target. 
This method results in positive identification of the signal and effective suppression of the 
background since, for example, the two-body K® 2 decay identifies itself by the unique 
momentum of the 7r° when viewed in the Kl rest frame. Once the Kl momentum is 
known, a large fraction of the K° L — > iftvv phase space is available for detection and all 
the major sources of background become manageable. 

Figure [7] shows the measured quantities in the experiment. The energies and times of 
the photons {E 1 , T 7 ) are measured by the calorimeter and the directions are measured by 
the preradiator. The preradiator also accomplishes particle identification: The first layer 
of scintillator identifies whether it is a charged or neutral particle, and dE/dx and the 
track pattern discriminates photons from Kl and neutrons. 

The Kl beam is collimated into a plane (slit beam). The decay vertex is measured 
by extrapolating the photon direction into the beam plane. The vertex can be obtained 
from one photon direction alone; the direction of the other photon, the time difference 
between the two photons, which corresponds to the path length difference between the 
two photons, and the 7r° mass work as extra constraints. These extra constraints are 
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Figure 7: Quantities measured in the experiment. 



effective in removing non-Gaussian tails, particularly in the angular measurements by the 
preradiator. The Kl energy is measured by time of flight using on RF bunched beam and 
direction is given by connecting the production target position and the decay vertex. 



2.2 Photon veto 

E787 has achieved a 7r° detection inefficiency of 10" 6 at photon energies of 20-225 MeV 
using a lead/scintillator detector. The inefficiency for the lowest energy photons is ~ 10~ 2 
and appears to be mainly limited by sampling fluctuations. The inefficiency for higher en- 
ergy photons is ~ 10~ 4 and appears to be limited by sampling fluctuations, shower escape 
and photonuclear reactions which may be contributing at comparable levels. Although 
the E787 group has attempted to study the origin of the residual inefficiency through 
measurements and simulations, considerable uncertainty remains. Thus, assuring the 
achievement of substantially higher photon detection efficiencies, for instance, at higher 
photon energies than the existing measurements, may be extremely difficult to establish 
reliably without actually performing the relevant measurements. The uncertainties are 
even larger for the region below 20 MeV photon energy which will be preferentially pop- 
ulated by K® 2 background events involving higher energy 7r°'s than observed in E787. 
It is for these reasons that, although advancements in detection efficiency may be possi- 
ble to achieve, we will only rely on small extrapolations from the E787 measurements in 
predicting the level of 7T° inefficiency. 

The goal in the K\ — > ir°uu experiment is to have 7r° detection inefficiency approach- 
ing 10~ 8 . This is feasible since both photons are generally in the higher energy range of the 
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E787 measurements where single photon detection inefficiencies of 10 -4 have been mea- 
sured. Since we also have kinematic handles available, we can suppress those kinematic 
configurations of K® 2 events with low energy missing photons and reasonably expect to 
achieve the goal. 

Because the momentum of the Kl is tagged, we can obtain the energy of the missing 
photons in K® 2 events by subtracting the measured energies of the two observed pho- 
tons from the Kl energy. Requiring significant total missing energy (i.e. (E(Kl) — 
E~f\ — E l2 ) as is generally the case for K\ — > 7T°vv events suppresses most poten- 
tial background events that contain lower energy missing photons (where the ineffi- 
ciency is greatest). However, in unusual cases when one of the missing photons has 
very high energy and one has very low energy an additional cut on missing mass (i.e. 

J (E(K L ) — E^i — E l2 ) 2 — (P(Kl) — P71 ~ P72) 2 ) is effective. Because the missing mass 
in K® 2 events is proportional to y/E miss i * E misS 2, where E miss is the energy of a missing 
photon, the missing mass also becomes small for the asymmetric energy sharing case. 
Figure [8] shows the missing mass vs. missing energy distribution of photons for K® 2 and 
K\ — > ifivv events. After removing the low missing mass and low missing energy region, 
we can suppress the low energy photons to achieve 10 -8 detection inefficiency for the two 
missing photons in K® 2 events. 
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Figure 8: Missing mass vs. missing energy distribution of photons for K® 2 odd, K® 2 even 
and K\ — > 7r°vi? events. 

The effect of small missing mass can be seen more directly by comparing E* distri- 
butions before and after the photon veto cut for K® 2 odd background events where one 
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photon from each ir° is missed (Figure 0). A dramatic peak above the E , * =230MeV after 
photon veto cut directly corresponds to the small missing mass []. This is one of the main 
reasons why the phase space below K® 2 peak (i?*o=249MeV) is used for the signal search 
region. 
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Figure 9: n° energy distribution in the Kl center of mass system (E* ) for K® 2 odd events 
before (top) and after (bottom) photon veto cuts. The vertical scale of the bottom figure 
is in the unit of 10 -8 events. 



3 Experimental design 



3.1 An overview 



Figure 10 shows a layout of the proposed experiment. A 500 fisr solid angle neutral 



beam is extracted at 45° to produce a "soft" Kl spectrum peaked at 0.65 GeV/c. Kaons 
in the range from about 0.5 GeV/c to 1.3 GeV/c will be used. The vertical acceptance of 
the beam (0.004 r) is kept much smaller than the horizontal acceptance (0.125 r) so that 
effective collimation can be used to severely limit beam halos. Downstream of the final 
beam collimator is a 3.5 m long decay region which is surrounded by the main detector. 
Approximately 16% of the kaons decay yielding a decay rate of about 25 MHz. The beam 
region is evacuated to a level of 10~ 7 Torr to suppress neutron induced n° production. 
The decay region is surrounded by an efficient Pb/scintillator photon veto detector. A 

m 2 +m 2 — m 2 

1 There is a one-to-one correspondence between E* and missing mass: E* = — — — 2if£_ p or 
small missing mass, E* is large. 
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Proton Beam 




Figure 10: Layout of the K\ — > 7r°vv experiment. 



Cerenkov detector ("catcher") covers the beam hole at the end of the beam line, 15 m 
downstream of the detector. 



3.2 The K L beam 



In recent years the AGS has achieved new records of intensity for synchrotrons. The 
present SEB peak extraction current is 7 x 10 13 protons/pulse. Coupled with a high cur- 
rent micro-bunched beam, good duty factor and extended availability during the RHIC 
era,0 the AGS is the ideal accelerator site for rare neutral kaon decay experiments em- 
ploying time-of-flight. The spectrum of neutral kaons produced at 45 degrees shown in 
Figure [U] (left) was obtained using a code based on a fit to CERN charged particle 
production data||. The spectrum shape is consistent with large angle charged kaon pro- 
duction data (BNL-E802) fitted by semi-phenomenological analytical methods|4]] ||. The 
Kl production yield is consistent within a factor of 2 between these two estimations. 

Neutrons from the production target could limit the beam intensity to be tolerated. 
In particular, neutrons above the pion production threshold can cause accidental back- 
grounds. The neutron momentum spectrum for 45° is shown in Figure [TTJ (right). At this 
large production angle, the neutron flux is largely suppressed above 800MeV/c, the pion 
production threshold of the n-p interaction. 

The micro-structure of the beam provides further suppression of the neutrons. Fig- 



2 RHIC is projected to operate for 30 to 40 weeks per year and requires injection from the AGS for 2 
hours/day. Thus, approximately 22 hours/day are available for AGS proton operation. 
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Figure 11: Momentum distribution of Kl and neutron fluxes 




Figure 12: Arrival time of K^s and neutrons with respect to photons at 10m downstream 
from the production target. 



12 



ure [12| shows the arrival time of K^s and neutrons with respect to photons at 10m from 
the production target. Neutrons with momenta between 0.9 and 2.4 GeV/c fall into our 
arrival time of interest (i.e. Kl with momenta between 0.5 and 1.3 GeV/c). Within this 
time window, the neutron to Kl ratio is improved by a factor of 5. Despite the fact that 
a low energy beam is used here, the effective n/K L ratio is as good as or better than in 
higher energy experiments. 



3.3 The detector 



In the forward detection region the primary photon detector system illustrated in Fig- 
ure O consists of two sections: a fine grained preradiator in which the photons are 
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Figure 13: K\ — > -n^vv detector. 



converted and the first e + /e~ pair is tracked, followed by an 18 radiation length (Xo) 
calorimeter in which the remaining energy of the photon shower is measured. The pre- 
radiator consists of 42 layers each with plastic scintillator, 0.03 X of heavy metal and 
dual coordinate drift chambers. The preradiator which has a total effective thickness of 
1.5 Xq functions to measure the photon positions and directions accurately in order to 
allow reconstruction of the Kl decay vertex while also contributing to the achievement of 
sufficient energy resolution. The calorimeter is constructed using thin Pb sheets formed to 
accept scintillating fibers in an arrangement similar to that recently made for the KLOE 
experiment at DA$NE; in our implementation of the calorimeter, a significantly higher 
visible light fraction will be used to further improve the energy resolution. In order to 
reduce the cost of the detector, we are also considering a shashlik calorimeter, which is 
made of layers of thin Pb sheet and plastic scintllator read out by WLS fibers running 
perpendicular to the layers. 

The barrel photon detector is made of layers of 1mm lead and 7mm extruded scintil- 
lator. The light from the scintillator is read out by embedded wave length shifter (WLS) 
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fibers. A prototype test with 4-m long WLS fibers shows a light yield of 12.1 p.e./MeV 
with a timing resolution of a =0.85nsec ||. Downstream of the main 7r° detector, a beam 
hole photon counter consists of layers of lead and lucite or aerogel Cerenkov detectors 
designed to be insensitive to neutrons. An initial prototype test using tagged photon and 
neutron beams shows promising results satisfying the experimental requirements 



4 Expected background levels and sensitivity 

The dominant background for K L —>■ ^vv is the CP violating decay Kl — > 7r°7r° (K n2 ) 
with a branching ratio of 9 x 10~ 4 . By tagging the Kl momentum as well as determining 
the energy and direction of 7s, one can reconstruct the kinematics. In the case of even 
pairing where one n° is missing, a kinematic cut on the monochromatic E* is effective. In 
the case of odd paring where one photon from each 7r° is missed, the 7r° mass requirement 
(m 77 ) is effective. An additional photon energy cut E* vs. \E* X — E* 2 \, where E* x and 
E* 2 are the energies of 7s in the Kl center of mass system, is also very effective in further 
suppressing the K® 2 background. 
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Figure 14: Distributions of E* vs. \E* X — E* 2 \ after tt mass requirement for the Kl — > 
7r 7r° (left) and for the signal (right). The solid line encloses the signal region. 

Figure [14] shows distributions of E* vs. \E* Y — E* 2 \ , after ir° mass requirement for the 
Kl — * 7T 7r° (left) and for the signal (right). The band at E* =249MeV corresponds to 
the even pairing background, which is suppressed by a E* cut at 225MeV/c as discussed 
above. The remaining band corresponds to the odd pairing background, which is confined 
into the band due to the constraint of 2 detected photons to form the n° mass. The solid 
line shows the signal region. As can be seen in the plots, we can obtain further background 
suppression at the cost of a modest acceptance loss. This redundancy in background 
rejection provides a safety margin in the experiment as well as a tool to positively identify 
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the signal. Table [I] gives the estimated acceptance factors for the K® 2 backgrounds. 



Requirement 




_0_0 „jj 

7r 7r oaa 


7r°7r° even 


7 combinations 


1 


4 


2 


solid angle 


0.36 


0.34 


0.37 




U. ( z 


U.Uo 


n 7/1 
U. 1 4 


TP* 


u.oo 


u. i u 


n nns 


photon veto 




1.3 x 10~ 8 


4.2 x 10~ 8 


Photon energy cuts 


0.35 


0.07 


0.17 


BR 




9 x 10" 4 


9 x 10" 4 


7 conv. and reconst. 


0.45 


0.45 


0.45 


Acceptance 


0.016 


3.5 x 10~ 14 


1.6 x 10" 14 



Table 1: Acceptance for K® 2 backgrounds. 



Other potential sources of background include neutron production of 7r°s, other Kl 
decays like K ei and Kl — > 77 and A — > mi decays. Suppression of most backgrounds 
is accomplished by the hermetic high efficiency photon detector along with kinematic 
constraints. A summary of the background estimates is given in Table [|. The signal 
is estimated to exceed the background by an order of magnitude with the backgrounds 
dominated by K® 2 . We expect that the actual background levels will be determined 
reliably from the data. 



Process 


Modes studied 


Main 


events 


K L -> n°vi> 






70 


K L decays (7) 


7r 7r ,7r 7r 7r ,7r°77 


7r°7r° 


7 


K L -> 77 


77 


77 


0.04 


Kl decays (chrg) 


71 ^ e T V , 71 ± jl T V , 7T + 7T _ 




0.01 


Kl decays (7, chrg) 


7r + 7T _ 7r ,7T ± / T ^7, 7r ± / =F ^7r ,7T + 7T~7 


7T + 7T~7r 


0.003 


Other decays 


A -> n°n, K~ -> tt'tt , S + -> n°p 


A — >■ 7r°ra 


0.03 


Interactions 


n, K L , 7 


n — > 7t° 


0.5 


Accidentals 


n, i^L, 7 


n, ATl, 7 


0.3 



Table 2: Estimated event levels for signal and background. 



Table |T| gives factors leading to the estimated acceptance of approximately 1.6%. It 
includes the solid angle, photon conversion and reconstruction factors and phase space 
(E* Q ) acceptance in addition to the photon energy cuts discussed above. The inefficiency 
due to accidental spoiling of good events is estimated to be < 10 % for a threshold of a 
few MeV and a timing window of 2 ns. The expected number of K\ — > tx vv events to 
be accumulated for 8000 hours of beam is ~ 50 assuming the SM central value for the 
branching ratio of B= 3 x 10" 11 . 
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5 Summary and prospects 



A detector system has been presented to measure the rare decay K\ — > -n^vv. The 
experiment is designed to allow definitive observation of a large sample of events with 
a signal that exceeds backgrounds by an order of magnitude. Special features of the 
AGS allow provision of an intense pulsed beam of neutral kaons suitable for time-of- 
flight measurements. The two independent signal selection criteria of photon veto and 
kinematics allow large background rejection and reliable measurement of the background 
levels. Significant redundancies and contingency factors are built into the technique. The 
standard model origin of CP violation will be confirmed and the complex phase parameter 
T) determined to a precision of < 15% or the absence of K\ — > ix^vv can be established 
at a level inconsistent with the standard model. 
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